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Nitric-oxide planar laser-induced fluorescence (NO PLIF) was used to visualize and measure 
centerline streamwise velocity of the Orion Crew Exploration Vehicle (CEV) Reaction Control 
System (RCS) Jets at NASA Langley Research Center’s 31-Inch Mach 10 Air wind tunnel. 
Fluorescence flow visualizations of pitch, roll, and yaw RCS jets were obtained using different 
plenum pressures and wind tunnel operating stagnation pressures. For two yaw RCS jet test 
cases, the PLIF visualizations were compared to computational flow imaging (CFI) images 
based on Langley Aerothermal Upwind Relaxation Algorithm (LAURA) computational fluid 
dynamics (CFD) simulations of the flowfield. For the same test cases, the streamwise velocity 
measurements were compared to CFD. The CFD solution, while showing some unphysical 
artifacts, generally agree with the experimental measurements. 
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0 2 = oxygen molecule 

P = pressure (Pa) 

P = reaction control system jet chamber pressure (Pa) 

P 0 = tunnel stagnation pressure (Pa) 

Q x = spectroscopic notation for transitions where parity of upper and lower states is positive and where AJ = 0 

Q = quenching rate (1/s) 

S f - = fluorescence signal 

T = temperature (K) 

T 0 = tunnel stagnation temperature (K) 

V L = velocity in laser direction (m/s) 

V s = velocity in streamwise direction (m/s) 

p = density (kg/m 3 ) 

(j) = fluorescence yield 

cj no = quenching cross-section of NO (cm 2 ) 

O' n 2 = quenching cross-section of N 2 (cm 2 ) 

o 0 2 = quenching cross-section of 0 2 (cm 2 ) 

t no = NO lifetime at tunnel stagnation temperature and stagnation pressure (s) 

X NO = molar concentration of NO 

Z Ni = molar concentration of N 2 

Xo 2 = molar concentration of 0 2 

Sv D = Doppler shift (1/cm) 

Sv L = laser detuning from line center (1/cm) 

Sv p = pressure, or collisional, shift (1/cm) 

Av d = Doppler broadening line width (1 /cm) 

Av l = laser spectral line width (1 /cm) 

Av p = pressure broadening linewidth (1 /cm) 

v 0 = centerline frequency (1 /cm) 

I. Introduction 

C OMPARATIVE experimental studies are necessary for the qualification and development of accurate 
computational fluid dynamics (CFD) models. CFD aids experimental vehicle configuration development during the 
iterative design process by lowering overall costs and reducing the turn-around time on designs. On the other hand, 
CFD models need to be validated against experimental data over a wide range of conditions to justify assumptions and 
verify model constants. This need is amplified when the simulations involve complicated physical processes that are 
typical of flows associated with hypersonic flight vehicles. The forebody flow of such vehicles include the effects of 
ablation, real gas effects, and thermal and chemical nonequillibrium, while the wake flow include flow separation, shear 
layers, and plumes. The transition to turbulence in both forebody and wake flows of these vehicles further complicates 
CFD modeling. In this paper, we consider such a wake flow: the aftbody flowfield of the Orion Crew Exploration 
Vehicle (CEV) upon reentry with the use of Reaction Control System (RCS) jets. 

In ground test facilities used to simulate reentry, few suitable spatially and temporally resolved flowfield 
measurement techniques exist. In this study, nitric oxide planar laser-induced fluorescence (NO PFIF) imaging and 
molecular-tagging velocimetry (MTV) are used to obtain three-dimensional visualizations and centerline streamwise 
velocities of the RCS Jets, respectively. The intensity of any given pixel in a PFIF image depends on many local flow 
parameters in addition to information about the optical system setup, complicating the extraction of any single local 
property. By modeling the fluorescence signal, direct comparisons between CFD predictions and the experiment can be 
made. This method is known as computational flow imaging (CFI) and has been applied to study similar hypersonic 
flows. 1 " 4 In this study, qualitative comparisons are made between CFI and experimental fluorescence images of the 
RCS jets interacting with the wake flow of the Orion CEV. The MTV velocity measurements with known uncertainties 
provide a quantitative means to compare the CFD and PFIF centerline streamwise velocities. The purpose of this paper 
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is both to report PLIF visualizations and velocity measurements of the Orion CEV RCS Jets and to evaluate a more 
detailed approach for comparing PLIF and CFD results. 


II. Experiment and Procedures 

A. Test Facility 

The experiments were performed in the 31 -Inch Mach 10 Air Tunnel facility at NASA Langley Research Center. 
The facility has a nominal Mach number of 9.7 with a nominal stagnation temperature, T {h of 1,000 K. Three large UV- 
grade fused silica windows, transparent in the ultraviolet down to approximately 190 nm, form three walls (top, side and 
bottom) of the 31 -inch square test section and a model injection system forms the fourth side wall. The top window 
provides access for the laser sheet to interrogate the flow, and the side window provides access to image the 
interrogated flow. Two facility stagnation pressures, Pq, were investigated: 350 psi (2.41 MPa) and 1,300 psi (9.30 
MPa). Further details of the flow properties at these conditions can be found in Hollis et al. 5 or from the authors. Run 
durations lasted for approximately 1 minute. For a thorough description of the facility, please see the paper by Micol. 6 

B. Test Model 

The experiments were conducted on a 0.025 
scale (5 -inch diameter) Orion Crew Exploration 
Vehicle (CEV) model, shown in Fig. 1. It had 
inserts containing 10-degree half angle conical 
Reaction Control System (RCS) jets with 0.7 mm 
(0.0275 in.) throat diameter and with an area ratio of 
13.4. The RCS consisted of yaw, pitch and roll 
nozzle jet pairs. In the present paper, only data from 
single-nozzles are presented. Jet pairs were also 
studied during the wind tunnel tests. The jet fluid 
composition was maintained at approximately 4% 

NO and 96% N 2 . Three model jet chamber 
pressures, P c , were investigated: -45 psi (0.31 
MPa), -132 psi (0.91 MPa), and -425 psi (2.93 
MPa). The axis of symmetry of the model was 
oriented at an 18° AoA with respect to the free 
stream. Reference 7 contains more information 
about the model. 

C. Laser System 

The laser system consisted of three main 
components: a pump laser, a dye laser, and a 
frequency-mixing unit. The second and third 
harmonics of a pulsed Spectra Physics Pro-230 
Nd:YAG laser, at wavelengths of 532 nm and 355 nm respectively, comprised the pump laser. The 532 nm output was 
used to pump a Sirah Cobra Dye laser to achieve a 622 nm output. The 622 nm output of the dye laser was then mixed 
with the 355 nm output of the pump laser in a Sirah Frequency Conversion Unit to produce -5 mJ per 10 ns pulse at 
226.256 nm. This wavelength was chosen to excite strongly fluorescing spectral lines of NO near the Q x band head. In 
this notation, Q denotes that the rotational quantum number, J, is the same for both the upper and lower states of the 
transition, i.e, A/=0. The subscript 1 indicates that the parity of both the upper and lower states involved in the transition 
is positive, meaning that in both states, the electronic spin is aligned with the total angular momentum of the molecule. 

To form the laser sheet, the collimated 226.256 nm output was passed through a 36 mm focal length cylindrical lens, 
oriented to diverge the beam in the direction parallel to the free stream. It was then passed through aim focal length 
spherical lens that collimated it in the freestream direction and focused it in the spanwise direction, forming a -75 mm 
wide by -0.5 mm thick sheet in the test section. For experiments involving molecular tagging velocimetry (MTV), 
multiple lines were required to tag NO. Therefore, a filter (a fine-toothed hair comb) was placed between the optics and 
the test section to form -1 mm line widths (FWHM) with -2 mm gaps. Further details of the laser system can be found 
in Bathel et al. 8 



Figure 1: Orion Command Module RCS Test Model (606D- 
TSP) assembled view with numbered pressure ports and jet 
locations (from Reference 7). 
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D. Camera System 

For the flow visualization experiments, a Princeton Instruments PI-MAX II with a gated image intensifier and a 
512x512 pixel CCD array was used. Custom filters (with <1 % transmission at 226 nm and >80 % transmission at 235- 
280 nm, made by Layertec GmbH) were used to transmit the fluorescence signal while rejecting the laser scatter. 
Images were acquired using a 45 mm F/1.8 UV lens (Cerco Model type 2073, manufactured by EADS Sodern), 
achieving a spatial resolution of about 2.7 pixels/mm. NO fluorescence was imaged at a viewing angle normal to the 
laser sheet at a rate of 10 Hz, matching the laser pulse rate, with flow-stopping time resolution (-1 ps). 

For the MTV experiments, a Cooke Dicam Pro with an intensified 1280x1024 pixel array interline progressive scan 
CCD was used. Only a 1280x288 pixel region was used so that an image read-out rate of 10 Hz, corresponding to the 
laser pulse frequency, could be achieved. A Nikon 105 mm F/4.5 UV-Nikkor lens was used, achieving a spatial 
resolution of 17.5 pixels/mm. The camera was operated in double shutter mode, wherein the camera was capable of 
acquiring an image pair with a 500 ns temporal separation (between the end of the first gate and the beginning of the 
second). The first intensifier gate duration was 20 ns, the minimum duration possible for this camera. Once a gas 
molecule has been tagged (excited by the laser), the fluorescence intensity decays exponentially in time, where the 
exponential time constant depends on the local pressure, temperature and gas composition. The Cooke Dicam Pro does 
not allow for independent gain settings for the first and second exposure. Therefore, in order to have a measurable 
second exposure, it was necessary to use the minimum gate delay (500 ns) and to keep the second gate open long 
enough (300 ns) to collect a sufficient number of photons. Furthermore, it was necessary to attenuate the first exposure 
by using the minimum gate duration setting (20 ns) and by advancing the camera gate 1 5 ns with respect to the laser 
excitation to create a shorter effective gate of 5 to 10 ns. These settings resulted in comparable intensities between the 
first and second exposures. The camera and the laser are triggered by a LabSmith LC880 Programmable Logic 
Controller to create a relative delay between the laser excitation and first exposure. Fine adjustments to the first gate 
duration were made by varying co-axial cable length between the trigger and the camera, yielding delay adjustment 
increments on the order of 1 ns (whereas the camera software allows delay adjustments only in coarse increments of 20 
ns). In summary, the effective duration of the first gate was between 5-10 ns, the delay between the end of the first gate 
and the beginning of the second was 500 ns, and the duration of the second gate was 300 ns. For a thorough description 
of the velocity measurement technique and more details about the camera setup, please see Bathel et al. 8 

E. PLIF Flow Visualization 

The -75 mm wide laser sheet is directed into the top tunnel window with the sheet oriented parallel to the freestream 
and normal to the camera system, which views the model and flow from the side window. The laser sheet excites the 
seeded NO molecules in the jet fluid. The resulting fluorescence is then imaged by the camera system. The laser sheet 
was scanned spatially through the flow using optics attached to a translation stage to visualize different cross-sections of 
the jet fluid. During such scans, the camera is fixed, so the magnification of the images varies by several percent. To 
account for the change in magnification, as well as lens and perspective distortion, a series of reference images were 
taken of a dotcard, which is a sheet marked with equally spaced square dots at a known position. 9 Each run investigated 
a single actuated RCS jet operated at two jet chamber pressures and a constant tunnel stagnation pressure. The test 
sequence is shown in Table 1. Further details on the PLIF flow visualization experimental procedure can be found in the 
paper by Danehy et al. 10 


Table 1. Test Matrix. Two conditions (denoted by an A or B appended to the Run number) were tested in each run. 


Run # 

Run Type 

RCS Jet 

P 0 (psi) 

Pc (psi) 

4A 

PLIF Visualization 

Pitch 

350 

45 

4B 

PLIF Visualization 

Pitch 

350 

136 

3A 

PLIF Visualization 

Pitch 

1300 

145 

3B 

PLIF Visualization 

Pitch 

1300 

425 

5A 

PLIF Visualization 

Roll 

350 

44 

5B 

PLIF Visualization 

Roll 

350 

130 

6A 

PLIF Visualization 

Roll 

1300 

144 

6B 

PLIF Visualization 

Roll 

1300 

420 

10A 

PLIF Visualization 

Yaw 

350 

45 

10B 

PLIF Visualization 

Yaw 

350 

132 

11A 

PLIF Visualization 

Yaw 

1300 

144 

1 IB 

PLIF Visualization 

Yaw 

1300 

424 

17A 

PLIF MTV 

Yaw 

350 

44 

17B 

PLIF MTV 

Yaw 

350 

131 
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F. PLIF Molecular Tagging Velocimetry 

In order to perform velocimetry, a slight modification is made to the laser sheet used for flow visualization: after 
optics form the beam into a sheet, a fine-toothed hair comb is placed between the optics and the test section to form ~1 
mm line widths with ~2 mm gaps, as described above in Laser System. The laser comb is directed vertically downward 
with the -75 mm side being parallel to free stream and normal to the camera system, as the laser sheet is in PLIF flow 
visualization. For the velocimetry data reported herein, the laser comb was not translated spatially but rather remained 
fixed at the jet center so repeated frames could be acquired. As in flow visualization, the laser excites the seeded NO 
molecules in the jet fluid. The resulting fluorescence is sequentially imaged twice by the camera system, where the first 
acquisition occurs at excitation and the second acquisition occurs after the excited lines of NO convect downstream for 
500 ns, as described above in Camera System. The single velocimetry run reported herein acquired multiple image pairs 
at the jet centerline. This run investigated the yaw RCS jet while operating at two jet chamber pressures and a single 
tunnel stagnation pressure. To account for lens and perspective distortions, reference dotcard images were again taken at 
several locations in the flowfield. For details of the run parameters, please see Table 1. Further details on the PLIF MTV 
experimental procedure can be found in the paper by Bathel et al. 8 

III. Data Analysis and Computational Methodology 

A. Flow Visualization Image Processing and Analysis 

Images are first processed to remove distortions associated with the imaging system. The UnwarpJ function of the 
National Institute of Health’s ImageJ program is used to apply this correction. An image acquired of the dotcard 
located at the model centerline is spatially transformed using UnwarpJ so that the dots are equally spaced and are 
aligned in vertical and horizontal rows to match the physical dotcard. The transformation is then applied to the dotcard 
images taken away from the model centerline and to the PLIF images to remove lens and perspective distortion. The 
background fluorescence is then subtracted from the PLIF images and a false color map is applied. In the Virtual 
Diagnostics Interface (ViDI), the transformed dotcard images are positioned relative to the Orion model according to 
their known offsets from the jet centerline. The dotcard images taken away the model centerline are scaled to match the 
dot size, placement and spacing of the actual physical dotcard. Linear functions are fit to the position and scaling 
information in such a way that the image offset from the jet centerline dictates relative position and scale. This 
procedure corrects for the variation in magnification caused by scanning the laser sheet position. Subsequently, PLIF 
images are imported into ViDI and aligned with the model with accuracy slightly less than 1 mm. The PLIF 
visualizations are positioned and scaled according to their recorded offsets and the linear functions obtained from the 
dotcard images, resulting in a three-dimensional representation of the jet. Further details on using ViDI to display PLIF 
images can be found in the paper by Inman et al. 11 

B. Molecular Tagging Velocimetry Image Processing and Analysis 

In brief, images acquired during the second camera gate are compared to the corresponding images acquired during 
the first camera gate. The displacement of the fluorescing lines of NO as a function of the relative time delay between 
the images is used to determine the flow velocity perpendicular to the lines. A 1 -dimensional cross-correlation method 
is used to determine the physical shift of the lines. For a thorough description of the MTV processing and 
corresponding error analysis methods used for the Orion CEV RCS Jets, please see Bathel et al. 8 

C. Computational Fluid Dynamics 

The CFD simulations shown in this paper were performed for a range of wind tunnel conditions and RCS nozzle 
flow conditions using Langley Aerothermal Upwind Relaxation Algorithm (LAURA) code. 12,13 Calculations used a 
laminar flow model with air as a perfect gas both in the free stream and in the nozzles. Nozzle grids extended to the 
throat, allowing some capture of the internal flow gradients. The nominal tunnel free stream conditions were used. The 
nozzle throat conditions were calculated assuming a 1 -D isentropic air nozzle flow. Calculations used a half-symmetric 
domain to save computational resources. This was thought adequate in the present case of the yaw nozzles, which are a 
sufficient distance from the centerline that the results should not be affected by the introduction of the artificial 
symmetry boundary. Grids contained from 6.5M to 52M points, though the solution used here was computed on a 13M- 
point grid. The standard LAURA bow shock and boundary layer adaptation was used, and the capsule’s wake grid was 
kept as short as is practical to keep up the grid density. Calculations used both local and global (time consistent) time 
integration with no noticeable difference in near-plume phenomena. 

For velocity comparisons, streamwise velocity was extracted from the CFD after a coordinate transformation. For 
the CFD computations shown in this paper, both the wind tunnel fluid and the RCS jet fluid were assumed to be air, 
which was treated as a perfect gas. For the NO PLIF CFI calculations shown in this paper, the entire fluid domain was 
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assumed to be composed of the experimental jet fluid ( Zo 2 = %n 2 = 0.96, % N o = 0.04). Ideally, the wind tunnel fluid 

and the RCS jet fluid could be modeled as two distinct gases whose compositions could be varied independently. 
However, the CFD computations were performed with a single fluid, air, for both the wind tunnel and RCS jet fluids. 
CFI computations required that the composition contain some NO in order to perform calculations of fluorescence 
intensity (i.e., the fluid could not be assumed to be air, as in the CFD computations); this suggested two possible 
prescriptions for the composition: 1) air with some seeded NO, or 2) the experimental jet fluid composition of N 2 
seeded with NO. Prescription 1 was tested; however, the results showed significant, unphysical amounts of intensity 
decay near the nozzle exit (a high-pressure area) due to 0 2 quenching. Since the experimental jet fluid did not contain 
0 2 (and since NO and 0 2 can react at high-pressures), a significant level of quenching only occurs in regions where the 
jet fluid has mixed with the wind tunnel fluid, which contains 0 2 , and not in the core of the jet plume. Thus, prescription 
1 was a poor approximation of the experiment. Prescription 2 was also tested, and although it fails to account for 0 2 
quenching in these mixing regions far from the jet centerline, the resulting jet-plume structures more closely matched 
the experimental visualizations; results presented in this paper therefore used this second prescription. The assumption 
of a single composition (N 2 seeded with NO) as both the freestream and RCS jet gas did complicate the CFI processing, 
as described in the following section. 


D. Computational Flow Imaging 

The intensity at any given pixel of a PLIF image has a functional dependence on the local temperature, pressure, 
velocity, density, and species, as well as the laser setup. Equation (1) gives a simplified expression for the fluorescence 
intensity’s dependence on the local state variables, where Sf is the fluorescence intensity, $ is the fluorescence yield, 
Xno is the local mole fraction of NO, p is the local gas density, f B is the Boltzmann fraction, B X2 is the Einstein 
absorption coefficient, g is the spectral overlap integral, and where the summation occurs over all probed transitions, i. 
There are approximately six physically significant transitions of NO near 226.256 nm, namely the strongly fluorescing 
spectral lines near the Q\ band head. For the calculations herein, these are the only transitions that are investigated. 
Equation (1) is relevant to flows where NO is the probed gas and where there is negligible saturation of the molecular 
transitions and rapid repopulation of the laser-excited ground state. 1,3,4 

Sf*<l>-XNO-p-Y}fB- B n-g\i ( 1 ) 


Equation (1) does not directly show how the fluorescence intensity depends on the local state parameters such as 
temperature and pressure since several of the terms are themselves functions of local state parameters, so the variables 
used within the relation will be further described herein. The quantity (f> is simply the ratio of the number of photons 
emitted as fluorescence to the number of photons absorbed, given by Eq. (2), where A 2 \ is the Einstein spontaneous 
emission coefficient and Q is the quenching rate. The quantity A 2 \ is given by 1 /t no , where t no is the NO lifetime at the 
tunnel stagnation temperature and pressure. 


. ^21 
^21 + 2 


( 2 ) 


The quantity Q in Eq. (2) is the sum of the quenching rates due to all molecular species present, as shown in Eq. (3), 
where k is the Boltzmann constant, m N0 is the molecular mass of NO, P is the local gas pressure, T is the local gas 
temperature, <j no is the NO quenching cross-section, Xo , i s the local mole fraction of 0 2 , m 0 is the molecular mass of 
0 2 , cr 02 is the 0 2 quenching cross-section, x Nl i s l° ca l mole fraction of N 2 , m N2 is the molecular mass of N 2 , and 
<j n is the quenching cross-section of N 2 . 3 


l 8-h-T 

\x-m NO 


P 
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m A 
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m A 


m n 


m A 


(3) 


The quenching cross-sections are temperature dependent; however, a theoretical representation of the dependency is 
unavailable, necessitating the usage of empirical models. The quenching cross-section models used within this study, 
shown in Eq. (4), were taken from Settersten et al. 14 

[34.8 -e 48 - 5/r + 96.0 -e-^j-lO- 20 (4) 

< 7 0 , = [22.5 -e 6h9,T +7.41- e 38° o/ r].i 0 2 ° 

<j Ni = [ 0 .OO 8 + 2.04 • e ~ 2250IT + 96.0 • e ~ uim,T ]. i o -20 


The quantity f B of Eq. (1) is given by Eq. (5), where h is the Plank constant, c is the speed light, B v is the effective 
rotational constant for the NO molecule which includes the effects of vibration/rotation interaction, and J is the 
rotational quantum number of the transition; development of this representation is found within the text by Eckbreth. 15 
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(5) 


h-c-B r n 

Equation (6) defines the quantity g, where g L is the spectral laser profile (Lorentzian or Gaussian profile, depending 
on laser), g a is the spectral absorption profile (Voigt profile), and v is the frequency. The Voigt profiles is a convolution 
of the Lorentzian-shaped pressure broadening profile with the Gaussian-shaped Maxwellian velocity distribution, which 
is associated with Doppler broadening. The spectral profile representations in Eq. (6) are developed in the text by 
Eckbreth. 15 Evaluation of Eq. (6) requires numeric integration and/or a lookup table. 


S = \g L 'ga dv 
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For all representations, Av D is the Doppler broadening width, Av P is the pressure broadening width, v L is the center laser 
frequency, Av L is the laser linewidth, v 0 is the center absorption frequency, dv L = v L - v 0 is the laser detuning from 
linecenter, dv P is the pressure shift and dv D is the Doppler shift. 

To reduce computation time for volumetric calculations, the quantity g a was approximated by assuming that the 
pressure broadening, usually assumed to have a Lorentzian profile was instead a Gaussian profile of the same width. 
This double-Gaussian simplification to the Voigt profile has an analytic solution. Furthermore, if the quantity g L was 
represented by a Gaussian profile as well (the so-called triple-Gaussian assumption), the additional convolution is 
analytically solved. The resulting analytic representation of the quantity g, shown in Eq. (7). 3 

f . ; „ w,' [Sv P +Sv D +Sv L f 


g(v,) = 


4 ■ ln(2) 


Av p+Av d +Av l 


7T-[ AVp +Avl + A v 2 l 


Expressions for Av D and 8v D are known exactly and are given by Eq. (8), where V L is the gas velocity parallel to the 


laser beam’s direction of propagation. 


Av d = 


2-v 0 1 2 • ln(2) -k-T 


s: V 0' V L 
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C 

The quantities Av L , and 8v L are determined from the experimental setup, whereas the quantities Av P and 8v P are 


modeled by Eq. (9), where the pressure P is in Pa. 
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The validity of triple-Gaussian spectral profile is explored in Experimental Results and Comparison with CFD , 
where the CFI fluorescence signal at the jet centerline is presented for the Gaussian- Voigt spectral overlap integral, the 
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spectral representation using the Gaussian laser profile and the Voigt absorption profile, at the highest jet pressure. Also 
compared in Experimental Results and Comparison with CFD is the simplified CFI expression shown in Eq. (10). This 
expression is a two dimensional polynomial fit for T and P and includes several assumptions. The assumptions include a 
rapid repopulation of the ground state, negligible saturation of the fluorescence transitions, a triple-Gaussian spectral 
overlap integral representation (assuming Zo 2 = 0? Zn, = 0.995, Zno = 0-005, v L = v 0 , and 8v£>= 0 for broadening and 
shift computations as well as quenching computations), a temperature range between 100 K and 500 K, a pressure range 
between 2 Pa and 30,000 Pa, and that Sf is linear with Zno-* A notable assumption is that quenching from 0 2 was 
neglected. At the time of formulation of Eq. (10), the updated quenching cross-sections of Eq. (4) were unavailable, so 
the quenching cross-sections from Paul et al. 17 were used. Further information on the simplified CFI representation can 
be found in Wilkes et al. 4 

S fmkes oc Xno -[-0.0018-r- 2 - 12 P 1 + 46.0-T~' 8 9 p\ (10) 

Please note that the remainder of this discussion will be focused on the CFI representation presented in Eq. (1), 
assuming the simplified g presented in Eq. (7). 

The CFD procedures assumed that both the jet and the ambient were a single fluid composition, which was air. In 
order to adapt the CFD solution to an acceptable form for NO PLIF CFI, all calculations assumed that the entire 
computational domain consisted of the exiting jet fluid ( Zo 2 = 0, Zn 2 = 9-96, Zno = 0.04). This assumption ignores 0 2 
quenching and results in a uniform computed background fluorescence outside the jet. 

All CFI calculations were done in Tecplot 360™ CFD post-processing and visualization software. In order to 
reconstruct a three-dimensional CFI jet, the computational domain was sliced into a series of vertical planes, oriented 
parallel to the tunnel streamwise direction, like the laser sheet. The planes were taken to extend past the jet plume in the 
spanwise direction to approximate the full jet plume structure with a distinct number of planes, a procedure akin to 
scanning the laser sheet through the jet plume for NO PLIF. The CFD results are used to determine the computational 
fluorescence signal for each of the sliced planes. Each plane then has the background fluorescence (typically just a few 
% of the maximum signal) subtracted off to account for the false assumption that there was also NO in the ambient wind 
tunnel fluid. While this subtraction does affect the ratios of fluorescence intensities elsewhere in the images, it did not 
significantly change fluorescence intensities in the jet fluid. The planes are then cropped to only show the jet plume 
region. This procedure was used to remove high levels of artificial fluorescence from the bow shock at the forebody of 
the Orion CEV that resulted from assuming NO was everywhere within the domain. The planes then have the same 
false color maps applied as the experimental PLIF images. Because the computed fluorescence intensity is numerically 
proportional to the actual intensity recorded by the intensified CCD, the CFI images are scaled so that the pixel intensity 
at the jet exit of the CFI centerline image matches that of the PLIF centerline image. The CFI images are then imported 
into the ViDI environment, where the planes are placed with respect to the jet centerline plane based on their 
coordinates obtained from Tecplot 360™. 


IV. Flow Field Structure 

Figure 2 shows several of the flow features of the Orion CEV flow field, which is typical of hypersonic capsule-type 
blunt body flows. Hypersonic blunt body flows are characterized by presence of a bow shock, which decelerates and 
compresses the oncoming flow. After passing through the shock wave, the gas expands outward in the radial direction 
towards the edge of the capsule’s heatshield, causing the local Mach number to increase. The pressure difference 
between the forebody and aftbody drives the flow to expand supersonically around the model’s shoulder (junction 
between the heat shield and the aftbody). Under some conditions the flow will separate and form a supersonic mixing 
layer, or shear layer, a boundary between the external flow and lower-speed recirculating flow. For example, on the 
leeward side of the model, the flow separates on the aftbody, as shown in Fig. 2. Finally the outer flow is processed by 
a recompression shock wave, turning the gas downstream. The size of the recirculation zone is sensitive to free-stream 
properties, model shape, vehicle AoA , gas chemistry, levels of turbulence and surface roughness, to name several (but 
not all) import factors. Flow on the windward side of the model can stay attached to the model surface, as shown in 
Fig. 2. This gas is also processed by a recompression shock downstream of the capsule. 
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Figure 2. Schematic of a capsule-type reentry vehicle in a hypersonic flow field (adapted from Ref. 13). Black 
lines indicate shock waves and streamlines as well as pressure iso-bars. The iso-bars show the area of expansion 
around the leeward and windward cone shoulders. 


In addition to the aforementioned flow structures around a blunt body shape, the presence of RCS jets on the 
vehicle’s aftbody adds another layer of complexity to the flow. Although the jets interact with both the wake flow and 
with each other, studying an isolated jet in crossflow can aid in the interpretation of the results presented later. Figure 3 
shows side and isometric views of underexpanded jets in supersonic crossflows, obtained from the literature. 18,19 A bow 
shock wave forms ahead of the jet due to the low streamwise momentum of the jet fluid relative to the crossflow. The 
interaction of the bow shock wave with the upstream boundary layer causes a the flow to separate forward of the bow 
shock. Downstream of the bow shock wave, spanwise vorticity in the boundary layer rolls up around the jet, forming a 
horseshoe vortex. At the lip of the jet orifice, Prandtl-Meyer expansion occurs due to the large jet-to-crossflow pressure 
ratio resulting in an underexpanded jet. The jet gas is subsequently compressed by a barrel shock wave and later by a 
Mach disk. A shear layer forms between the barrel shock wave and the bow shock wave, which provides the mechanism 
for the production of turbulent eddies. A counter-rotating vortex pair is shown in the isometric view of Fig. 3. The 
vortex pair interacts with the boundary layer and convects high momentum, high temperature gas from the freestream 
down to the wall surface. 




Figure 3. Schematics of an underexpanded supersonic jet in crossflow. The side-view (a) is an instantaneous 
snapshot of the flow field (Ref. 18). The isometric view (b) is a time-averaged sketch of the flow field (Ref. 19). 
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V. Experimental Results and Comparison with CFD 
A. PLIF Flow Visualization 

Figures 4 and 5 show a summary of the NO PLIF flow visualization data set obtained in the wind tunnel tests. They 
compare pitch, roll and yaw jet flows for low tunnel stagnation pressure (Fig. 4) and high tunnel stagnation pressure 
(Fig. 5) cases. While these results are qualitatively similar to those described in detail in Ref. 11, these images were 
obtained with a slightly different Orion model configuration in a different test entry. (In the current paper, the Orion 
configuration is 606D and the test number was Test 455). The main focus of this paper is the comparison between the 
CFD computations and the fluorescence imaging and velocimetry for the yaw jet case. The experimental results in Figs. 
4 and 5 will therefore only be briefly described. 



Figure 4. Isometric views of experimental PLIF volumetric flow visualizations with P 0 = 350 psi. Column (a) has 
P c ~ 45 psi. Column (b) has P c ~ 132 psi. Figures (al) and (bl) show the pitch jet; (a2) and (b2) show the roll jet; 
(a3) and (b3) show the yaw jet (note model rotation 90° about the sting in the third row relative to the first two 
rows). 
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Figure 5. Isometric views of PLIF volumetric flow visualizations with P 0 = 1,300 psi. Column (a) has P c ~ 132 psi. 
Column (b) has P c ~ 425 psi. Figures (al) and (bl) show the pitch jet; (a2) and (b2) show the roll jet; (a3) and 
(b3) show the yaw jet (note model rotation). 


In the top row of Fig. 4, the pitch-up jets issue into the centerline wake flow region, which is a large separated 
region having relatively low static pressure. Near the jet exit, there exists a large separated flow region, containing 
slowly moving gas, resulting in very weak crossflow. Strong fluorescence signal continues along the jet’s trajectory 
until the jet reaches the shear layer between the separated wake flow region and the flow passing over the heat shield 
and over the shoulder of the model. The jet fluid is then observed to flow mostly downstream, though some jet fluid 
impinges on the aft-body thermal protection system as well. The roll jets, in the middle row of Fig. 4, initially 
propagate vertically (nearly azimuthally) out of the jet exit. When the jet fluid reaches the oncoming crossflow passing 
around the side of the model, the jets are deflected downstream towards the sting. The roll jets appear to partially form 
the counter-rotating vortices shown in Figs. 2 and 3. To study the yaw jet, the model was rotated 90° counter clockwise 
(as viewed from upstream) about the sting. The yaw jet visualizations shown in the bottom rows of Figs. 4 and 5 also 
have been rotated by 90° compared to the pitch and roll visualizations. Considering cases with similar jet chamber 
pressures, P c , the pitch jets appear largest in diameter and the yaw jets appear smallest. These trends are caused by the 
spatial pressure distribution on the backshell as well as the varying size of the separated flow region at the different jet 
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locations. Higher pressures and smaller separated regions reduce the extent of the jet flow. The higher P c in the right- 
column figures results in a larger diameter jet, evident in pitch, roll and yaw jets. 

In Fig. 5, the stagnation pressure, P 0 , increases by a factor of 3.7 while the values of P c increase by about a factor of 
3 from the corresponding cases in Fig. 4. Thus the ratio of P c /P 0 stays roughly constant in Figs. 4 and 5. The main 
factor that changes is the Reynolds number, with Fig. 5 having a 3.7 times higher Reynolds number than Fig. 4. In 
comparing these two cases, significant changes in the flow features (size, shape, trajectory of the jets) are not observed. 
The fluorescence intensity decreases more quickly with downstream distance, for example in the roll-jet cases, but this 
is probably because the higher pressure results in faster quenching of the NO fluorescence and consequently lower 
signal. This does not imply any fluid mechanical significance. One would expect more turbulent flow features in the 
higher Reynolds number cases, but the lower Reynolds number case already appears rather turbulent at the edges of the 
jets observed with this type of volumetric visualization. Numerous other views including single-plane images and 
animations of these flow visualizations have been produced and are available from the authors. 

B. PLIF Velocity Measurements 

Figure 6 shows the average experimental velocities and measurement uncertainties for two chamber pressures, P c , 
obtained during the same run. The lower chamber pressure case generally has larger uncertainties than the higher 
chamber pressure case due to a smaller number of measurements averaged together. 8 Generally, the uncertainties in the 
velocity increase with downstream distance and increase dramatically near the edge of the jet. The larger uncertainties 
near the boundary of the jet are due primarily to the unsteady nature of the flow features. Velocity data near the jet exit 
were removed due to a systematic error in the velocity calculation caused by a large, negative quenching gradient over 
the width of the line. 8 



Figure 6. PLIF MTV yaw jet measurements with P 0 = 350 psi. The position of the center of the black line relative 
to its local green scale bar indicates its mean velocity. The thickness of the black line relative to its local scale bar 
indicates the uncertainty in the velocity measurement. Figure (a) has P c ~ 45 psi; (b) has P c ~ 132 psi. 


C. CFD Temperature, Density, Pressure, and Velocity 

The CFD temperature, density, pressure, and streamwise velocity maps for yaw jet operation shown in Fig. 7 are all 
taken in the plane passing through the jet centerline position and oriented in the vertical plane, parallel to the streamwise 
direction. The temperature images in Fig. 7 show several of the flow features described in Figs. 2 and 3. The smooth 
change in temperature caused by the expansion fan from the shoulder of the model is observed in the top left corner of 
each image. This high speed gas then interacts with the jet plume, resulting in an oblique bow shock, observed as a 
diagonal high-temperature region in Figs. 6(a) and (b). Some unphysical oscillations are observed in this bow shock, as 
well as in the downstream jet plume. These artifacts may be caused by non-optimal grid density and alignment in the 
wake region. 13 The post-shock gas cools as it mixes with the cold RCS jet fluid, which has a stagnation temperature of 
approximately room temperature and has been expanded past Mach 4. A room- temperature stagnant region exists 
towards the bottom-right side of both Figs. 6(a) and (b). The computed density, pressure and velocity presented in Fig. 
6 corroborate with this description and are qualitatively similar to many (time-averaged) aspects of Figs. 2 and 3. 
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Figure 7. CFD temperature, density, pressure, and streamwise velocity results with P 0 = 350 psi. Column (a) has 
P c ~ 45 psi. Column (b) has P c ~ 132 psi. 
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D. CFI and PLIF Flow Visualization Comparison 

Figures 8 and 9 compare theoretical (CFI) and experimental PLIF fluorescence images for the yaw Orion RCS jet. 



Figure 8. CFI and Experimental PLIF yaw jet flow visualizations with P 0 = 350 psi and P c ~ 45 psi. Column (a) 
shows the CFI images. Column (b) shows the comparable PLIF flow visualization, (al) and (bl) are side views; 
(a2) and (b2) are bottom views; (a3) and (b3) are back views; (a4) and (b4) are isometric views. 
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Figure 9. CFI and Experimental PLIF yaw jet flow visualizations with P 0 = 350 psi and P c ~ 132 psi. Column (a) 
shows the CFI images. Column (b) shows the comparable PLIF flow visualizations, (al) and (bl) are side views. 
(a2) and (b2) are bottom views. (a3) and (b3) are back views. (a4) and (b4) are isometric views. 

The comparison in Fig. 8 and 9 is made by simultaneously displaying images from many closely spaced vertical slices 
through the yaw jet. The images show that the CFI is able to convey very similar three-dimensional geometries of the 
yaw RCS jet for both of the chamber pressures. The ‘back’ views in the third row of both figures show the expected 
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counter-rotating vortices forming on the edge of the jet as in a classical jet in a crossflow, shown in Fig. 3(b). These are 
highlighted by arrows in Fig. 9(a3) and (b3). The darker appearance of the CFI PLIF in the bottom two rows of each 
figure is not fluid-mechanically significant. It is an artifact of the color table used to display the images and the spacing 
of the CFI-PLIF planes. An obvious difference between the CFI and experimental plumes is that the CFI plumes are 
smooth, owing to the laminar nature of the computation, while the experimental data show turbulent structures on the 
edges of the jets (though they also show laminar flow in the core of the jets, as shown in Fig. 10). The origin of these 
turbulent structures on the outer edges of the jet is described above, in the explanation of Fig. 3(a). 



Figure 10. CFI and Experimental PLIF yaw jet flow visualizations with P 0 = 350 psi. Column (a) shows the CFI 
images. Column (b) shows the comparable PLIF flow visualizations, (al) and (bl) are bottom views with P c ~ 45 
psi. (a2) and (b2) are bottom views with P c ~ 132 psi. 

Whereas Figs. 8 and 9 show the three-dimensional structure of the yaw jet by displaying data from multiple planes 
simultaneously, Fig. 10 shows some of the finer structures of the jet as revealed by single image planes. The CFI and 
PLIF centerline images, shown in Fig. 10, illustrate some discrepancies between theory and experiment. The CFI has a 
more diffuse and smaller sized barrel-shock boundary. The furthest downstream extent of the jet, where one might 
expect a Mach disk to appear, shows a rounded shape in the CFI while the experimental image shows a pointed 
structure, indicated by the arrow in Fig. 10(b2). The presence of barrel shock waves and Mach disk boundaries are 
expected in the PLIF images since sudden changes in pressure and temperature are closely related to changes in 
fluorescence (Eq. 10). Some small, unphysical diagonal features are observed inside the barrel shock in Fig. 10(a2). 
These may be grid alignment artifacts, 13 suggesting that additional refinements to the grid might help improve the 
comparison between theory and experiment. 
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E. CFI Representation Comparison 

Figure 11 shows the results of a comparative study conducted on the center plane of the yaw jet flow visualization 
test case with P 0 = 350 psi and P c ~ 132 psi. Figure 1 1 compares two different NO seeding concentrations as well as 
two different ways to compute and display the CFI PLIF images. Figure 1 1 (a 1 ) was computed assuming that the 
spectral overlap integral, g, is Gaussian- Voigt (a Gaussian laser profile and a Voigt absorption profile) and the mole 
fraction of nitric oxide was Xno = 0.04; Fig. 1 1 (b 1 ) assumes g is triple-Gaussian and Xvo = 0.04; Fig. 1 l(a2) assumes g 
is triple-Gaussian and Xvo = 0.005; and Fig. 11 (b2) uses the Sf mikes representation. Figure (a3) is a line plot taken 
approximately down the centerline of the jet, downstream of the jet exit, showing relative intensities for the CFI 
representations in Figs. 1 l(al), 1 l(bl), 1 l(a2), and 1 l(b2). The fluorescence profiles in Fig. 1 l(a3) are scaled relative to 
each other to approximately match the fluorescence intensity downstream. Also shown in Fig. 1 l(a3) is a cutoff for the 
Sf wnkes representation used to saturate the jet exit and increase the surrounding intensity, as shown Fig. 1 l(b3). Figures 
1 1 (al) and 1 1 (b 1 ) look nearly identical, demonstrating that the triple-Gaussian spectral representation approximated the 
Gaussian- Voigt profile very well. This is a significant result because it means that the convolution integral in Eq. 6 does 
not have to be computed. This reduces the computational time by 2-3 orders of magnitude (depending on the number of 
integration increments used). This comparison is also shown more clearly in Fig. 1 l(a3) where a single profile through 
the center of the jet is shown. 



Figure 11. CFI yaw jet flow visualizations and line plot with P 0 = 350 psi and P c ~ 132 psi. (al), (bl), (a2), (b2) 
and (b3) are CFI images of the same fluid domain using different assumptions, explained in the text. (a3) is a line 
graph comparing how the different models indicate the same profile down the center of the jet. 
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Figures 11 (a2) and (b2), computed with nearly a factor of 10 lower NO seeding, though identical overall jet mass 
flow rate, show a more diffuse and slightly smaller-appearing jet boundary than the higher seeding cases in Fig 1 l(al) 
and (bl). The apparent decrease in jet size is probably due to the color table, which is scaled to the maximum value in 
the image: higher peak values near the jet exit make the jet appear narrower (the CFD computation was the same for 
both figures). Figures 1 l(a2) and (b2) show that the simple Sf mikes expression computes a PLIF distribution (b2) that is 
very close to the more sophisticated triple-Gaussian computation in (a2). The fluorescence profiles shown in Fig. 
11 (a3) more clearly illustrate the effect of using a smaller concentration of NO. At lower concentrations, the 
fluorescence intensity increases dramatically near the jet exit, where there is the highest pressure and highest 
concentration of NO. At higher NO concentrations and higher pressures, the self-quenching of NO reduces the intensity 
near the jet exit. A drawback of using the lower NO concentration is that there is a very large change in signal between 
the jet exit and the downstream plume. So, if this concentration is used, the dynamic range of fluorescence intensity 
may be difficult to display on a linear color table as shown above (the jet exit will appear very bright and the 
downstream plume will not be visible). As a workaround, the Sf wakes intensity can be purposefully saturated (in the 
color table) near the jet exit, for example at the relative fluorescence intensity shown by the cutoff in Fig. 11 (a3), the 
results of which are shown in Fig. Il(b3). The CFI image presented in Fig. Il(b3) better shows the barrel shock 
boundary (as indicated by the white arrow) compared to 1 l(b2). The Sf wakes representation is much easier to implement 
and substantially faster than the representation of Eq. (1) assuming g is triple-Gaussian, with time savings on the order 
of 2-3 magnitudes. Also, it would be a simple matter to compute a new version of the Sf wakes expression for the higher 
seeding levels if desired. For future comparisons between computational and experimental PLIF images in RCS jet 
experiments, we recommend use of the simple Sf wakes expression. While it does not have all the physics embedded in 
the more complicated models, it is much easier to use, much faster, and can easily be implemented by a non-expert in 
Tecplot. Also, the level of discrepancies observed between the more sophisticated model and the experimental data are 
larger than the differences between Sf wakes and the more sophisticated models, at least for the RCS jets considered in the 
present study. 

F. CFD and MTV Velocity Comparison 

Experimental MTV measurements are shown overlaid on maps of CFD velocity computations for two different jet 
chamber pressures in Fig. 12(al) and 12(bl). For each run condition, three measurement locations have been chosen for 
a more detailed comparison. In Figs. 12(a2)-12(a4) and 12(b2)-12(b4), lines showing the CFD results are plotted along 
with MTV data, including measurement uncertainty error bars, for these selected locations. As Fig. 12 illustrates, the 
CFD velocities generally showed similar growth and decay trends as the MTV measurements, with the theory generally 
fitting within the window of uncertainty of the measurement. Near the edges of the profile, the disparity between MTV 
and CFD velocities increased significantly. Also, the disparity between theory and experiment grows with distance 
downstream, as shown generally in Fig. 1 1 (b 1 ) . In the future, the larger uncertainty bars could be reduced by a factor of 
2-3 by averaging 4-9 times more measurements. This number of measurements (up to 900) could be obtained in a 
single tunnel run using the existing instrument. However, some uncertainties in the experiment 8 will not be reduced by 
averaging. 
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Figure 12. CFD and MTV centerline streamwise velocity comparisons with P 0 = 350 psi. Column (a) has P c ~ 45 
psi. Column (b) has P c ~ 132 psi. (al) and (bl) display the velocities using a color table, where the lines indicate 
the MTV velocity measurements and the background indicates the CFD velocity measurements. The arrows 
within (al) and (bl) indicate which profiles are displayed in (a2)-(a4) and (b2)-(b4), respectively. Within the 
profiles, the bars extending from the MTV data indicate the measurement uncertainty. 
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VI. Discussion and Conclusion 

This work summarizes experimental and computational PLIF flow visualizations of the Orion CEV RCS jets. It 
compares PLIF MTV velocity measurements with CFD solutions at the same conditions. Furthermore, a more detailed 
model for CFI PLIF was presented and compared with a simpler model used in past work. 

The CFI and PLIF flow visualizations generally agreed in the jet trajectory and shape. The disparities between the 
theoretical and experimental results can be attributed partly to the CFD solution and partly to the assumptions inherent 
in the CFI representation of Eq. (1). There could also be inadequate modeling of the NO laser-induced fluorescence 
process. A likely candidate for the differences between computations and experiment is the assumption that there was 
no 0 2 quenching in the flowfield; this could not be avoided in the present work because the CFD solution assumed that 
air was both the test gas and the RCS jet fluid. Had a unique jet fluid been used in the CFD solution, the model 
presented herein would be able to account for the effects of 0 2 quenching. Other possibilities for the discrepancies 
include saturation of the molecular transitions from high laser intensities, laser absorption by NO, and the time 
dependence of the fluorescence process, all of which were neglected in the current work. A simplified method of 
computing the PLIF intensity from CFD, given by S/ t mikes was shown to work satisfactorily and is recommended for use 
in the future. This curve-fitting approach could be extended to include, for example, different concentrations of NO, 
and possibly 0 2 . Such a model can be distributed to the CFD community to do their own comparisons quickly without 
knowledge of the underlying spectroscopy. The streamwise velocities of the CFD were mostly within uncertainty of 
our MTV measurement technique. While not providing a formal validation of the CFD codes, these comparisons give 
confidence that the main aspects of the RCS jet computation (jet size, shape, trajectory and core velocity) are being 
adequately modeled. 
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